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Abstract 



The function Q(x) := ]^n>i(V n ) sin(x/n) was introduced by Hardy and Littlewood [10] in their 
study of Lambert summability, and since then it has attracted attention of many researchers. In 
particular, this function has made a surprising appearance in the recent disproof by Alzer, Berg and 
Koumandos [1] of a conjecture by Clark and Ismail [2]. More precisely, Alzer et. al. have shown 
that the Clark and Ismail conjecture is true if and only if Q(x) > —ir/2 for all x > 0. It is known 
that Q(x) is unbounded in the domain x £ (0,oo) from above and below, which disproves the Clark 
and Ismail conjecture, and at the same time raises a natural question of whether we can exhibit at 
least one point x for which Q(x) < —tt/2. This turns out to be a surprisingly hard problem, which 
leads to an interesting and non-trivial question of how to approximate Q(x) for very large values of 
x. In this paper we continue the work started by Gautschi in [7] and develop several approximations 
to Q(x) for large values of x. We use these approximations to find an explicit value of x for which 
Q(x) < -vr/2. 
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1 Introduction 



Our main object of interest is the function Q(x), defined as 

«(.):= (5). (i) 

n>l 

It is easy to see that the above series converges for all x G C and that Q(x) is an odd entire function of 
exponential type one. This function was named "Flett's function" by van de Lune [15] and Doring [5], 
however later it was given the name "Hardy-Littlewood" function by Alzer et. al. [1] and Gautschi [7], 
and we will follow the latter convention in our paper. 

The function Q(x) was originally introduced in 1936 by Hardy and Littlewood [10], who have used 
it to construct a certain counter-example in their investigations of Lambert summability. In particular, 
Hardy and Littlewood have established that Q(x) = fLf_(y lnln(x)), which is just a short notation for 

liminf Q(x) / \/\ia ln(rr) > 0. 

Hardy and Littlewood also note that the "O-problems" for Q(x) are much like the corresponding problems 
for the Riemann zeta function £(s) on the vertical line Re(s) = 1 (note that it is known that |£(1 +it)\ = 
f2 + (lnln(t)), see Theorem 8.5 in [14]). 

It seems that this similarity in the behavior of Q(x) and ((1 + it) is the main reason why Q(x) 
has attracted so much attention after the work of Hardy and Littlewood. Flett [6] has proved that 

3 1 

Q(x) = 0(ln(x)4(lnln(x))2 +<E ) and that the same estimate is true for £(1 + it). Segal [13] has discovered 
many interesting properties of Q(x), such as the following identity 

J Q{X)AX = 25> = f - \ + £ (^)» J, ( 2v ^) , (2) 

q n>l rra>l 

where J\(z) is the Bessel function of order one, see [8]. As an application of (2), Segal proves that the 
Cesaro, Abel and Borel means of Q(x) are all equal tt/2. 

Codeca [3] has studied oscillation and almost periodicity properties of Q(x) and some other related 
functions. He also makes an important observation that Q(x) is very similar to some functions which 
appear as error terms in the estimation of the mean of certain number-theoretic functions. For example, 
if we define a "divisor function" cr_i(n) := ^2 d \ n l/d and S-i(x) := J2 n<x o'-i(n), then it is known (see 
[16], p. 100) that 

S^x) = ^-1 \n(x) -Y^U Q + 0(1), (3) 

n<x 

where (j)(x) := {x} — | and {x} denotes the fractional part of x. Note that both sin(a;) and (fr(x) are 
periodic functions, which means that the two functions Q(x) and S(x) := Y2 n <x(^/ n ) < i ) ( x / n ) should have 
many similar properties. This turns out to be a beneficial way of looking at these functions: among 
many other results in [3], Codeca has proved that both Q(x) and E(x) are B 2 almost periodic function 
and are unbounded from above and below. Codeca also mentiones that Delange [4] has extended the 
results by Hardy and Littlewood [10] and by Flett [6] and has proved that 

Q(x) = Q ± (VhMx)), Q(x) = 0({ln{x))%), x^+oo. (4) 
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Again, the latter result should be compared with the best known bounds £(x) = 0((ln(a;))3) (see [16], p. 

2 

88) and C(l + it) = 0((ln(t))5) (see formula 6.19.2 in [14]). The study initiated by Codeca was continued 
by Petermann [12], in particular he gave another proof of (4). 

Van de Lune [15] has performed the first numerical study of the function Q(x), in particular he has 
calculated many real and complex zeros of this function. Another method for computing zeros of Q(x) 
was developed by Doring [5]. 

Recently, the Hardy-Littlewood function Q(x) has appeared rather unexpectedly in the dispoof of 
the Clark and Ismail conjecture by Alzer, Berg and Koumandos [1]. Clark and Ismail [2] have studied 
functions $ m (x) := — x m il)^ m \x), where ip(x) := T'(x)/r(x) is the digamma function. They have proved 
that (x) is completely monotone on (0, oo) for 1 < m < 16, and they conjecture that this should 
be true for all m > 1. However, this conjecture was disproved by Alzer, Berg and Koumandos [1] by 
showing that it is true if and only if Q(x) > — | for all x, the latter statement being false in view of (4). 

Finally, we would like to mention the paper [7] by Gautschi, which was the main inspiration for 
our current work. Gautschi has developed two algorithms for computing Q(x) numerically: the first 
algorithm is based on the summation by quadrature and the second (more efficient) algorithm is based 
on truncating the series in (1) at n = [x\ (the integer part of x) and approximating the tail of this series. 
It is clear that this algorithm requires 0(x) arithmetic operations in order to obtain a single value of 
Q(x), thus it becomes impractical if x is very large. 

Our main results in this paper are several new asymptotic approximations for Q(x). Our first result 
gives an approximation for Q(x), which is extremely accurate in the domain x > 2000 and which requires 
only 0(x2 +e ) arithmetic operations. The second approximation is somewhat less accurate, but it requires 
only 0(x3) arithmetic operations. As an application of these two approximations we find that 

Q(8203872394818031742687.4 x tt) = -1.5970415 . . . , (5) 

which provides an explicit example to the result by Alzer et. al. [1] and answers the question raised by 
Gautschi in [7]. 

This paper is organized as follows. In Section 2 we review the approximation developed by Gautschi 
[7] and present our main results, Theorems 1, 2 and 3. In Section 3 we perform several numerical 
experiments: we investigate the accuracy of our approximations, we study the extremes of the function 
Q(x) and discuss the computations that led to the discovery of (5). The proofs of all results are presented 
in Section 4. 

2 Main results 

Let us introduce the notations that will be used throughout this paper. Given two functions / : D i— > R 
and g : D \-t (0, oo), we will write f(x) = 0(g(x)) (or equivalently, f(x) <C g(x)), if there exists an 
absolute constant C > such that \f(x)\ < Cg(x) for all x G D. When the constant C is not absolute 
but depends on parameters Pi, ■ ■ ■ ,p n , we will write f(x) = Plv .. iPn (g(x)) or f(x) •C Pll ... )Pn g{x). We will 
also use the notation f(x) ~ g(x), which stands for f(x) -C g(x) and g(x) <C f{x). Bernoulli numbers 
will be denoted by B n . Finally, for iGl, [^J := max{n < x : nfZ} will denote the integer part of x 
and {x} = x — [x\ will denote the fractional part of x. 

Our first result is an algorithm which allows to compute Q(x) to arbitrary precision in 0{x) arithmetic 
operations. This algorithm is a simple generalization of the approximation developed by Gautschi (see 
Section 3 in [7]). 
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Proposition 1. (0(x) algorithm) For any integer N we have 



N-l 



n=l fc>0 V y 

where c&(iV) are defined as follows 

Ck (N) := N 2k+1 n~ 2k ~ 2 - (7) 
For large N the coefficients c^(iV) can 6e computed via the asymptotic expansion 



m=l 



The proof of the Proposition 1 is rather simple, we just sketch the main steps and leave all the 
details to the reader. In order to obtain (6), we expand each term of the tail J2 n >N^/ n ) sin(x/n) in the 
Taylor-Maclaurin series and interchange the order of summation. Asymptotic expansion (8) follows at 
once by applying the Euler-Maclaurin formula to (7). 

As we will see in Section 3, for large values of x a good choice of N in Proposition 1 is N = |_^J . 
Therefore, the above algorithm needs 0(x) arithmetic operations and it becomes impractical for very 
large values of x. Our first main result is an asymptotic approximation to Q(x), which requires only 
0(xz +e ) arithmetic operations and is extremely accurate in the domain x > 2000. In order to present 
this algorithm, we will need to define the sine integral function 

Si(s) = / ™<&6t. (9) 



See Section 8.23 in [8] for many properties of this function. We will need only one of these properties, 
namely that for large values of x the sine integral can be computed via the asymptotic expansion 



M /„,si M 



Si(x) = - - cos(x) £(-l) fe |*| - sin ^ B- 1 )"^^ + Om(x-™~% M E {-1, 0, 1, . . . }.(10) 

k=0 k=0 

The above asymptotic expansion can be easily derived from formulas 8.215 and 8.233.1 in [8]. 

Theorem 1. (0(x^ +e ) algorithm) Assume that < e < |. Define G(x) := | sin and N = N(x) := 
[x^ +t \. Then 

N-l 

«*> = £;;' fa (s) + 23f' ta (£) + Si (£) + < u > 

n=l 



where £n(x) has the following asymptotic expansion: for M £ {0, 1, 2, ... } 

M 

X 



R f(2m-l) (N\ 

= - E(ggr T J' J + o M „^'"^% (12) 
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The proof of Theorem 1 will be presented in Section 4. This result should be compared with the 
approximation to the Riemann zeta function ((s) in the critical strip < Re(s) < 1, see Theorem 4.11 
in [14]. 

Note that in its simplest form the Theorem 1 gives us the following approximation: for < e < | 

Q(x) = £ - sin (|)+| + °^" l+£ )- (13) 

n=l 

To obtain (13) one should take M = — l{M = 0}in the asymptotic expansion (10) { resp. (12) } and 
combine these results with the formula (11). The constant tt/2 in (13) is reminiscent of the result by 
Segal [13] (which was already mentioned in the Introduction), namely that the Cesaro, Abel and Borel 
means of Q(x) are all equal tt/2. Formula (13) raises the following two natural questions: 

(i) Would it still be correct in the limiting case e = 0? 

(ii) Can we reduce the number of terms in the sum in the right-hand side of (13) to \_x a \ with some 
a<\! 

Our next result provides the answers to both of these questions. 

Theorem 2. (0(xa) algorithm) Assume that M and N are real numbers, such that M > 1, N > 1 and 
2ttMN 2 = x. Then 

i 

^= E n Sin (D + i + E (^) 3 smQ + 2v / 2^)+0(iV- 1 ln(x)+x^ivl). ( U ) 

l<n<iV l<m<M 

The proof of Theorem 2 will be presented in Section 4. This result should be compared with the 
approximate functional equation for the Riemann zeta function, see Theorem 4.13 in [14]. 

Note that if M > 1 is a fixed constant (which does not depend on x), then the second sum in the 
right-hand side of (14) is 0(x~i), and the number of terms in the first sum is N = LvV(2ttM)J. This 
provides the answer to question (i): we can take e = and iV w yfx in (13), but then we can only 
be certain that the error is 0{x~± ln(x)), and not 0(x~a) as (13) would suggest. Another important 
observation is that the parameters M and N are linked through the condition 2tcMN 2 = x, and if we 
decrease the number of terms in the first sum in the right-hand side of (14) we will at the same time 
increase the number of terms in the second sum. It is clear that we obtain the best order of approximation 
and the smallest number of terms in both sums in (14) if we take M = N = ^/x/(2vr), in which case the 
error becomes 0(x~^ ln(a;)). In particular, this gives the answer to question (ii). 

Remark 1. Segal [13] has raised the following interesting question. If we formally differentiate the 
identity (2), we would obtain 

nmx). (15) 

m>l 

The problem is that we do not know whether the series in the right-hand side of (15) converges. The 
connection between formulas (14) and (15) is provided by the following asymptotic expansion (see formula 
8.451.1 in [8]) 

7iJ (2V2nmx) = (^~) * sin (f + 2 ^ 27imx ) + (H^) » ( 16 ) 
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which implies that (14) could serve as a possible interpretation of (15). 

As we have discussed above, the choice M = N = y / x/(2ir) in Theorem 2 is optimal in the sence 
that it gives the highest order of approximation and the smallest number of terms in the two sums in the 
right-hand side of (14). If we take N ~ x a then the condition 2nMN 2 = x would imply that M « a; 1-2 ", 
thus is seems that it would be impossible to reduce the number of terms in (13) to [x a \ with some a < |. 
Surprisingly, this is not the case. Our next result states that the number of terms in (13) can be reduced 
to [x a \ with any a < |. 

Theorem 3. (0(x e ) algorithm) Let < e < | and define S = 5(e) := e2~e . Then 

q(x)= J2 ^ in © + i +a(ar5) - (17) 

l<n<x f - 

The proof of Theorem 3 will be presented in Section 4. To the best of our knowledge, there is no 
analogous result for the Riemann zeta function (unlike the previous two Theorems). 

3 Numerical results 

Our first goal is to investigate the accuracy of the approximations provided by Theorems 1 and 2. We will 
use Proposition 1 in order to compute the benchmark values of Q(x). We choose N = max(\_x\, 1000) 
and compute the coefficients Cj-(N) for < k < 9 using the asymptotic formula (8) with M — 6. Then 
we truncate the second series in the right-hand side of (6) at k = k max = 9. Experimenting with higher 
values of N, M and fc max we find that the above parameters guarantee the accuracy of at least 15 digits. 
The code was written in Fortran90 and we have used quadruple precision for all computations. 

We denote by Qi(x) the approximation to Q(x) obtained by setting e = 0.05 and M = 4 in Theorem 
1, and by Q2(x) the approximation to Q(x) provided by taking the optimal choice of parameters M = 
N = ^/x/(2ir) in Theorem 2. We present the results of our computations on figure 1. On figure la we 
see that Qi(x) provides an excellent approximation to Q(x): in the domain x > 2000 the absolute error 
is smaller than 10 -10 and in the domain x > 4500 the absolute error is smaller than 10~ 12 . In order to 
investigate the accuracy of the second approximation we would like to consider much larger values of x. 
Given the fact that for x large Qi(x) is very close to Q(x) and is much easier to evaluate numerically, 
we will use the values of Qi(x) as a benchmark. The results are presented on figure lb. We find that 
the error Q2{x) — Qi(x) is of the order 10~ 3 when x ~ 10 9 and of the order 10~ 4 when x ~ 10 12 . Note 
that these results are in perfect agreement with the error estimate 0(x~s ln(a;)) in (14). 

As our next goal, we have tried to find a value of x where Q(x) < —n/2. This would provide 
an explicit example for the key step in the disproof of the Clark and Ismail conjecture by Alzer et. 
al. [1]. Our first attempt was to look at the very large local maximums/minimums of Q(x). As was 
noted by Alzer et. al. [1], for every n e {0, 1,2, ... } the function Q(x) has a local maximum in the 
interval [2n7f, (2n + 1)tt] and a local minimum in the interval [(2n + 1)tt, (2n + 2)ir]. By checking all 
consecutive local maximums/minimums of Q(x) we have found several values of x, where Q(x) has a 
local extremum and the value of this local maximum/minimum is greater /smaller than the value of all 
other local maximums/minimums in the interval [0, x). The results of our computations are presented in 
Table 1. The lowest value of Q(x) that we were able to find using this approach is approximately —1.31. 
Given the fact that the function Q(x) grows very slowly (see (4)), it became clear that this brute force 
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(a) log 10 \Q(x) - Qi{x)\ for 50 < x < 5000 (b) Q 2 (x + x) - Qi(x + x) for < x < 5 x 10 4 

Figure 1: The error of the approximation. Here Q(x) is the value computed using (6) with iV = 
min(x, 1000) and M = 10; Qi(x) is computed using (11) and (12) with e = 0.05 and M — 4, and Q2^x) 
using (14) with M = N = y/x/(2ir). On figure (b) the black curve corresponds to x® = 10 9 and the red 
curve to xq = 10 12 . 



search method will not work. Therefore, we will have to bring in some other ideas in order to reach the 
value of Q(x) which is less than —n/2. 

In order to find these new ideas, we have looked at the proof of the estimate Q(x) = Q±(yln ln(x)) 
given in [1] and [10]. This proof is based on a rather interesting argument, which is both constructive 
and non-constructive. Let us summarize the main steps of this argument. We define a set of integer 
numbers Ai, such that q G M. if and only if all divisors of q satisfy q = l(mod 4). The first few elements 
of M. are 

M = {5, 13, 17, 25, 29, 37, 41, 53, 61, 65, 73, . . . }. 
Now, for each k > 1 we define an integer number 

K = K(k) := q, 

l<g<4fc+l 
q&M 

and for each j e {1,2,..., K} we define 

xf := (4j + xf := (4j + 3)k|. (18) 

In [1] and [10] it was established that there exist constants a G R, a G R and b > 0, b > 0, such that for 
every k there exist i = i(k) and j = j(k) such that 

Q(xf } ) > a + b^Mk), Q(xf ] ) < o - fev 7 !^)- (19) 

Using this result coupled with the trivial estimate K(k) < (4k + l) 4fc+1 it is not hard to establish that 
Q(x) = tt±(^/ln\n(x)) as x — > +oo. 

The intuition behind the definition (18) is rather simple: one can see that sia(x ( - /n) = 1 and 
sin^^/n) = —1 for all n which divide K. Since has many small divisors n, this shows that there 
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x/n 


12127812.568 


37324872.600 


50774112.608 


176438112.602 


413620512.598 


Q(x) 


4.3300 


4.3426 


4.3931 


4.4596 


4.4638 


x/tc 


3596987.431 


3841175.419 


51836087.411 


196661495.417 


580973087.402 


Q(x) 


-1.0512 


-1.1635 


-1.2406 


-1.2740 


-1.3134 



Table 1: Some "extreme" local maximums/minimums of Q(x) 



k 


K{k) 


max Q(x^) 

Ki<K 1 


max achieved at 


min Q(x^) 


min achieved at 


4 


1105 


4.1352 


i = 876 


-0.9262 


j = H03 


6 


27625 


4.2606 


i = 24966 


-1.1647 


j = 6785 


7 


801125 


4.4127 


i = 259627 


-1.2498 


j = 468148 


9 


29641625 


4.5752 


i = 7030990 


-1.4347 


j = 615949 


10 


1215306625 


4.6586 


i = 96330096 


-1.4717 


j = 69224831 



Table 2: Computing maximum/minimum values of Q(x\ ) and Q(xj ) 



will be quite many terms in the sum (1) defining Q(x\ ), which are equal to 1/n (or which are equal to 
—1/n, in the case of Q(x^)). Of course, it still takes a lot of work to deduce (19), as we have to show 
that the sum over all n which do not divide K is not too large. See [1] and [10] for all the details. 

Formulas (18) and (19) give us an algorithm for finding very large positive/negative values of Q(x): 

(k) ~(k) 

these extreme values will happen at points x\ and Xj . This is the constructive side of this result. At 

the same time, there is no information on how to find the indices i and j for which Q{x^) or Q{x^) 
would achieve the maximum/minimum values, therefore this part would still have to be done by brute 
force search, by checking all the indices in the range 1 < i < K(k). Clearly, this is only feasible if k is 
not too large, as K{k) grows very fast as k increases. 

The results of our computations are presented in Table 2. For k = 4 and k = 6 we have used the 
approximation Qi(x), provided by Theorem 1 with e = 0.025 and M = 4. For larger values of k we 
have used the approximation Q2(x), obtained from Theorem 2 by setting M = N = ^/x/(2tt). In 
this case the value of Q(x) at the extreme point was confirmed by computing it via the more accurate 
approximation Qi(x). All computations were performed on a regular desktop computer with Intel i7 
2600 quad-core processor, running Ubuntu Linux. In order to fully utilize all four cores of the processor, 
we have parallelized the Fortran90 code using OpenMP API. 

It is instructive to compare the results presented in Tables 1 and 2. First of all, one can check that 
the points where Q(x) attains extremely large local maximums/minimums in Table 1 are not located 
near the points x\ k ' or x^ defined by (18). Second, we see that the large values of Q(x) in Table 1 
happen at smaller values of x than those in Table 2. This provides a compelling numerical evidence 
that the result Q(x) = ^(^/ln ln(a;)) (which is derived via (19)) is suboptimal. A possible way to 
prove a stronger result would be to understand how to predict the points in Table 1, where Q(x) attains 
its extreme local maximums/minimums. Given the similarity between Q(x), ((1 + it), £(x) and other 
functions arising in Number Theory, this understanding could potentially lead to establishing stronger 
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3 
2.5 



1.5 



0.5 


-0.5 



2 
-2 

I* - 40 x* - 20 x* x* + 20 x* + 40 

Figure 2: A region where Q(z) < -§. Here x* = 8203872394818031742687.5 x vr. 

results about oscillation of all these functions. 

As we see from Table 2, there are no values of Q(xj) for k < 10 which are smaller than — tt/2. 
And k = 10 was the largest value that we could possibly check on our desktop computer in a reasonable 
amount of time: in order to perform the computations for the next value k — 13 (for which we need 
to evaluate Q(x) '), 1 < j < K(13) = 64411251125) we would need to wait for nine months. Luckily 
enough, after just two days of computations at the level k = 13, we have found that at the index 
f = 63683535496 we have 

Q( x *) = -1.5398404..., (20) 

where 

x * ■= x ^ = (4f + 3)AT(13) J = 8203872394818031742687.5 x vr. 
J 2 

The value in (20) is just slightly above — n/2. However, looking at the results in Table 1, we see that the 
local minimums of Q(x) tend to occur near points x such that x/tt = 0.4 (mod 1). We check the value 
of Q(x) at this point, and find that 

Q(x* -O.lvr) = -1.5970415.... (21) 

The value of Q(x* — O.ln) in (21) was computed using the second approximation Q2(x), provided 
by Theorem 2. Given the error estimate 0(x~s ln(x)) in Theorem 2 and our previous discussion of the 
numerical results shown in Figure 1, we would expect that Q2{x) — Q{x) ~ 10~ 7 when x ~ 10 22 . In 
order to confirm this and to check the accuracy of computation in (21), we have evaluated Qi(x* — O.ln), 
which is the approximation given by Theorem 1 with e = 0.025 and M = 4. We found that the difference 
between the two results is less than 10 -7 , as expected. Finally, in order to eliminate the possibility 
that there was some loss of accuracy in (21) due to the addition of a large number of terms in (14), we 
have evaluated Q2{x* — O.Itt) with the working precision of 60 decimal digits. We have used D. Bailey's 
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MPFUN multi-precision library for Fortran90. Again, the results were in perfect agreement with (21). 
While we cannot prove (21) rigorously (since we do not know the explicit values for the implied constants 
in (14)), these numerical results provide compelling evidence that (21) is correct, and therefore, provides 
an explicit example for the key step of the disproof of the Clark and Ismail conjecture by Alzer et. al. 

[1]- 

The plot of Q(x) near the point x* is presented on Figure (2). 

In conclusion, we would like to discuss the computation time of Q(x* — 0.l7r). This will help us to 
put the three approximations to Q(x) into perspective. It takes just five seconds to establish (21) using 
the 0(xz) algorithm provided by Theorem 2. Computing the same value using the 0{x^) algorithm 
presented in Theorem 1 took almost fourteen hours. By extrapolation, we find that the 0(x) algorithm 
from Proposition 1 would need more than forty million years in order to complete the same task. 



4 Proofs 



Proof of Theorem 1: We set h := 1/x and apply the Euler-Maclaurin summation formula to function 
G(x) = - sin f-J: for every j > 



h 



1 00 j 

1 -G(Nh)+ °( nh ) = G{u)du-^-^h 2k G^{Nh) (22) 

n>N+l J £ h k=l ^ >' 



+ 



h 2 i- 



(2J + 1)! 



B 2j+1 ({u/h})G^ +1 \u)du, 



Nh 



where B n (-) denote Bernoulli polynomials and B n = B n (0) are Bernoulli numbers. 

Changing the variable of integration u = 1/v and using (9) we find that the first integral in the 
right-hand side of (22) is equal to Si(x/N). 

Let us now estimate the second integral in the right-hand side of (22). Using induction it is easy to 
prove that there exists a sequence of polynomials {p n (u)} n >o, such that deg(p n ) = n and for each j > 
we have 



G& +1 \u) = u 23 - 2 \cos{u x ) Vw {u- x ) + sin^- 1 )^- 1 )] . 



(23) 



Formula (23) implies that for every j > there exists a constant C\(j) > such that for all u > we 
have 



\G^ +1 \u)\ < C.U) (a - " ■ a 2 ) 



(24) 



The periodic function |B2j+i({M//i})| is bounded from above by some positive constant C^j), therefore 
we obtain the estimate 



B 23+1 {{u/h})G^ +1 \u)du 



Nh 



< CiU) x C 2 (j) x / (u^~ 3 + u- 2 ) du = Oj ((Nh)-^- 2 ) . 



Nh 



Using the fact that Nh w x 2 +e as x — > +oo, we conclude 



h 2 ^ 
(2j + l)! 



B 2H1 ({u/h})G^ +1 \u)du = Oj (x~ {4j+2)e ) . 



(25) 



Nh 
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Combining (11), (22) and (25) we see that for every j > 

3 



S N (x) 



B 2k G^ 2k V( N x ) ( T -(*j+2)e\ 
2^(2k)\ x 2k + A 



(26) 



For every e G (0, |) and M > we can find j large enough so that (4j + 2)e > | + (AM + 3)e. Also, from 
(24) we obtain 



G (2*-l) 



2A- 



x- 2fc O fc ((AT/a;)- 4fc+1 ) = O fc (arH«-i)«). 



Combining these two facts with (26) we conclude that (12) holds. □ 

The proof of the Theorem 2 is based on the following two Lemmas. 
Lemma 1. (First and Second derivative test, see Lemmas 5.1.2 and 5.1.3 in [11]) 

Let f : (a, 0) i— > R be twice continuously differentiable and g : (a, (5) i— y [0, oo) be monotone. For 
i E {1,2} we define Aj = inf{|/^(-u)| : a < u < /3} and assume that A2 > 0. Then 



g{u) sin(/(u))dii 



< 16 x max{g(a+), g(/3— )} x min 



Ai a/As 



(27) 



Lemma 2. (Truncated Poisson summation formula, see Lemma 5.4.3 in [11] or Lemma 4.10 in [14]) 
Let f : [a, b] h->- R /icrae a continuous and decreasing derivative f'(u). Let g : [a, b] (->■ [0, 00) fre decreasing 
function, such that \g'(u)\ is also decreasing. Define a = f'(b) and (3 = /'(a). Then 



5> 

a<n<f> 



n e 



27ri/(n) 



E 



/ ( M ) e 2-i(/(«)-m U ) du + _ a + 2 )) + 0(|«/( a )|) (2? 



Q-|<m</3+i a 



Proof of Theorem 2: The proof will be carried in three steps. Our first goal is to establish the following 
identity: 

+ 0(iV~ 1 ln(x)). (29) 

Note that the conditions 2tcMN 2 = x, M >1, N >1 imply that 

M = 0(x), N = 0(y/x), x > 2tt. (30) 

We define a := N, b := x 2 , f(u) := — and g(u) := - and check that all conditions of Lemma 2 
are satisfied. We find that a = f'(b) = an d = f'{ a ) = M and note that for x > 2tt we have 
— 1 < a — ^ < 0. Applying Lemma 2 and taking imaginary part of both sides of (28) we obtain 

+ 0(N~ 1 ln(M)). (31) 



E- 



n>N 



X 

sin ' 

n \n 



E 



Q<m<M+\ 



CO 

/— sin ( — h 2nmu ) du 
u \u J 



N 



} — sin — 
^— ' n \n 

N<n<x 2 



E 



0<m<A/+| 



/ — sin ( — h 2nmu\ du 
u \u J 



N 
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Next, assume that m > 1 and w > x 2 . We define f(u) :— - + 2^mu and note that 

Ai = inf { I /'(w) I : x 2 < u < w} = ^ + 2irm, 

x 6 

A 2 = inf{|/"(u)| : x 2 <u<w} = \>0. 

x 5 

We apply Lemma 1 (the first derivative test) with the above function f(u) and g(u) := u^ 1 and conclude 
that for m > 1, x > 27r and any w > x 2 



/ — sin ( — h 27tmu ) du <C 
u \u / 
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Therefore for all x > 27r and m > 1 we have 



oo 

/ — sin ( — h 2nmu ) du — 
u Vu / x^ 



When m = we have the following trivial estimate 



1 /X 
— sin — 

It \U 



du < 



sm 



x 



du < 



x 1 
— du = -. 

U X 



The above two estimates and the fact that M = 0(x) give us 



E 

0<m<M+l \_ x 2 



— sin (— + 2i\mu ) du 
u \u 



M 1 1 

< -5- + - < 
x z x a; 



(32) 



Combining (30), (31), (32) and the following simple estimate 



E- 



X 

sm ' 

n vn 



«E 

n>x 2 



X 1 

77/ x 



we obtain (29). 

The second step consists in establishing the following result 



E 

0<m<A/+i 



N 



/ — sin ( — h 2i:mu\ du 



< N ln(x) + x 2 j\T2. 



(33) 



We change the variable of integration u — l/v and obtain 



N oo 

/ — sin ( — h 2nmu ] du = / - sin ( xu + 
u \u J J v \ 

AT" 1 



27rm 



du. 



(34) 
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Define f(v) := xv + 2Ttm/v and assume that < m < M, where M := [M — . Then for v > N 1 we 
have 



... , 2nm 9 
f (v) = x — > x - 2irmN . 



v 2 



Note that m < M < M - \ implies 



x - 2nmN 2 > x - 2ti { M - - \ N 2 



A) 2N 2 ' 

thus Ai = inf{|/'(f)| : A^ 1 < v < 00} = x — 2nmN 2 and Ai > 0. Applying Lemma 1 (the first 
derivative test) to the integral in the right-hand side of (34) we conclude that 

N 

16N 



f 1 . fx n \ . f 1 . / 27rm\ , 

/ — sin — h 2Txmu ) du = / - sin xv H dv <C 

J u \u J J V \ V J 



x — 2-KmN 2 



N 

Using the above estimate we obtain 

M 



m=0 



N 

— sin ( — h 27rmu ) du 



™ 16N " 16N 



< V ; -r- — ^ = V 7^ r < iV~ 1 ln(2 + M) < AT- 1 ln(a;), 

- x - 2vr(M - \)N 2 + 2tcjN 2 j^7iN(l + 4j) v ; v ; ' 

where we have changed the index of summation j = M — m and have used (30) and the fact that 

In order to finish the proof of (33) we have to consider the possible integral term with M < m < M+j 
in the sum in the left-hand side of (33). We again change the variable of integration v = 1/u and consider 
the following two integrals 

N 2N- 1 00 

f 1 fx \ [If 2vrm\ f 1 / 27rm\ . . 

/ — sin — h 2ixmu ) du = / - sin I xv H } dv + / - sin I xv H dv (36) 

J u \u J J v V v ) J v \ v J 

N- 1 2N- 1 

Define f(u) := xu + ^p, then f'(u) = x- 2 ^ and f"(u) = It is easy to see that 

A 2 = min{|/"(«)| : A^ 1 <u< 2N' 1 } = ^nmN 3 . 
Applying Lemma 1 (the second derivative test) we obtain the following estimate 

2N" 1 

/l / 27rm \ , , „ T . 1 1 1 , , 

- sin I ct H ] dw < (mN) 2 < x~*Ni (37) 

AT-l 

where in the last step we have used the facts that m > M — | and MiV = x/{2nN). Next, we find that 
for u > 2N- 1 

f'( M ) = x =- > x - -ir(M + -)A^ 2 = > -, 

J K J u 2 v A ! 4 8 2' 
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thus by Lemma 1 (the first derivative test) we obtain 



CO 



1 / litrn 

sin xu H ) du < Nx~ l < iV -1 , (3£ 



u \ u 



2N~ 



where in the last step we have used the upper bound Nx 1 = N 1 /(2ixM) <C N 1 . Equations (36), (37) 
and (38) show that for M < m < M + \ 



N 



/- sin (- + 2Tcmu] du <^ N 1 + x ^N*. 



The above estimate combined with (35) complete the proof of (33). 

Now we have all the ingredients for the last step of the proof of Theorem 2. We combine the two 
results (29) and (33) and conclude that 



n ( n ) 



n \n 

n>N 0<m<M+± 



oo 

/ — sin ( — h 2irmu ) du 
u \u J 



+ 0(N- 1 ln(x) +x- 1 *N^) (39) 



4 L 

^ + tt ^ J (2V27rmx) + 0(N' 1 ln(x) + x~*N? 



l<m<M-i 

where we have used the following integral identities (see formulas 3.721.1 and 3.868.1 in [8]): 



oc 



/ — sin ( —] du = — , [ — sin ( — h 2Tvmu] du = n Jo(2V2wmx) . 
u \uJ 2 J u \u J 



o 



In order to obtain (14) one has to combine (39), the asymptotic expression (16) for J (2^2 irmx) and 
the following simple estimate 



l<m<M+\ 



□ 



The proof of Theorem 3 is based on the following result from the theory of exponential sums. 

Theorem 4. (Theorem 2.9 in [9]) Let q > be an integer and Q = 2 q . Suppose that f has q + 2 
continuous derivatives on an interval I and that I C (N,2N]. Assume also that there is some constant 
F such that for l<r<q + 2we have 

\f( r \u)\ « FN~ r , uel. (40) 

Then 

e 2«/W < p^N 1 '^ + F^N. (41) 
The implied constant in (41) depends only upon the implied constants in (40). 
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Proof of Theorem 3: Let iV be any positive number. For u > N we define 

'X s 



S(u):= sin (j 



N<n<u 



and S(u) := for u < N. We also define f(u) := x/(2im) and F := x/N. Computing f^ r \u) one can 
easily check that for iV < u < 2N and 1 < r < q + 2 

r! x FN- r < \f {r) (u)\ < — x FN~ r , 



therefore all the conditions of Theorem 4 are satisfied. Taking the imaginary part of (41) we conclude 
that for all g > 0, iV > 1 and N <u < 2N 

1 i q-\-3 

S{u) < 9 ii^iY 1 -!^ + x -ijy 2 (42) 
where Q = 2 q . We apply the integration by parts and find that for all N < M < 2N we have 



n \n/ / it M + \ j u 



< N~ l max <^\S(u) | : iV - ^ < u < M + i j < g x 4 «- 2 JV 4Q + - 2 + x" 1 ^, 

where in the last step we have used (42). 
Let us define N := [x e \ and 

j = j{x) := [log^N-'x^l , a = a(q, e) := ((q + 3)e - 1)/(4Q - 2). 

Dividing the interval N < n < x 3//5 into subintervals of the form Ni < n < N 2 < 2Ni we obtain 

j_l 2 *+ 1 jV-l 

E-sin(-) = -sinf-) + —sin ( — J (44) 



V<n<x 3 / 5 »=0 n=2*iV 2iN<n<x 3 / 5 



< 9 x*e- 2 ^ 4 « + -2 (j + 1) + x _1 (x 3/5 - AT) < 9 a;"" ln(x) + af §, 

where we have used the facts that j <C ln(x) and A" rs x e . We choose g = [ e_1 J — 1 an d find that 

Le- 1 Je + 2e-l e i 

a = ^ — in,-, > -i — rrr-; >f2 < = 5(e). 

2L^ 1 J+i - 2 2L 6-1 J+ 1 - 2 w 

This fact combined with the estimate (44) shows that 

V -sm(-) =OJx- s ). (45) 



n \n 

x e <n<x 3 / 5 



Taking e = 1/10 in (13) we find that 

Q(?)= E ^ sin © + i +0(ar§) - 

l<n<x s / 5 

Combining this result with (45) ends the proof of Theorem 3. □ 
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